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Algebraic Fermi liquid from phase fluctuations: "topological" fermions, vortex 
"berryons" and QED 3 theory of cuprate superconductors 

M. FranzQ and Z. Tcsanoviejj] 

Institute for Theoretical Physics, 
University of California, Santa Barbara, CA 93106 
(February 1, 2008) 

Within the phase fluctuation model for the pseudogap state of cuprate superconductors we identify 
a novel statistical "Berry phase" interaction between the nodal quasiparticles and fluctuating vortex- 
antivortex excitations. The effective action describing this model assumes the form of an anisotropic 
Euclidean quantum electrodynamics in (2+1) dimensions (QED3) and naturally generates non-Fermi 
liquid behavior for its fermionic excitations. The doping axis in the x-T phase diagram emerges as 
a quantum critical line which regulates the low energy fermiology. 



Perhaps the most intriguing property of high tempera- 
ture superconductors is the anomalous character of their 
normal state [jij. This "strange metal" stands in stark 
contrast to the relatively benign features of the supercon- 
ducting phase which can be understood rather accurately 
within the framework of a d-wetve BCS-like phenomenol- 
ogy with well defined quasiparticle excitations 0. 

In this Letter we propose a theory of the pseudogap 
phase in cuprate superconductors based on the follow- 
ing premise: a successful phenomenology of the strange 
metal should be built by starting from a comprehensive 
understanding of the adjacent superconducting state and 
its excitations. The spirit of our approach is the tradi- 
tional one ||, |j| but turned upside down. Usually, the 
strategy is to first understand the normal state before 
we can understand the superconductor. In the cuprates, 
however, it is the superconducting state that appears 
"conventional" and its quasiparticles "less correlated" 
and better defined. Having adopted this "inverted" 
paradigm, we proceed to study the interactions of the 
quasiparticles with the collective modes of the system, 
i.e. fluctuating (anti) vortices (our strategy here is sim- 
ilar to that of Ref. g). We show that in d-wave super- 
conductors these interactions take a form of a gauge the- 
ory which shares considerable similarity with the quan- 
tum electrodynamics in (2+l)-dimensions (QED3). In 
the superconducting state, where vortices are bound, the 
gauge fields of the theory are massive and the low energy 
quasiparticles remain well defined excitations. This is the 
mundane Fermi liquid state in our inverted paradigm. In 
the normal state, however, as vortices unbind, our QED3 
like theory enters its massless phase and it abandons this 
"inverted Fermi liquid" protectorate in favor of a weakly 
destabilized Fermi liquid characterized by a power law 
singularity in the fermion propagator which we call alge- 
braic Fermi liquid. We compute the spectral properties of 
fermions in our theory and find that it captures some key 
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FIG. 1: Phase diagram of a cuprate superconductor. 



qualitative aspects of the available experimental data. 

We concentrate on the portion of the pseudogap phase 
above the shaded region and below T* in Fig. 1. We 
assume that Cooper pairs are formed at or somewhat be- 
low T* but the long-range phase coherence sets in only 
at the superconducting transition temperature T sc <C T* 
Q. Between T sc and T* the phase order is destroyed by 
unbound vortex-antivortex excitations of the Cooper pair 
field §@, §. Inthispseudo gap regime the d-wave super- 
conducting gap is still relatively intact j|, ||] and the dom- 
inant interactions are those of nodal quasiparticles with 
fluctuating vortices. There are two components of this 
interaction: first, vortex fluctuations produce variations 
in superfluid velocity which cause Doppler shift in quasi- 
particle energies if) . This effect is classical and already 
much studied [{jij, 01. Second, there is a purely quan- 
tum "statistical" interaction, tied to a geometric "Berry 
phase" effect that winds the phase of a quasiparticle as it 
encircles a vortex jl2], [l3] . It is this quantum mechanical 
interaction that ultimately causes the destruction of the 
Fermi liquid in the pseudogap phase. 

Our starting point is the partition function 

Z = J V^\y,t) J £>*(r,r) J Vip(r, t) exp [-S], 
S = J dr J d 2 r{^d T ^ + ^H^> + (1/ 5 )A*A}, (1) 
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where r is the imaginary time, r = (x, y), g is an effective 
coupling constant, and ^ — (V'f,V'j.) are the standard 
Grassmann variables. The Hamiltonian H is given by: 
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with 7i e = ^(p - f A) 2 - e F , p = -iV (we take 
Ti = 1), and A the d-wave pairing operator Jl2| , 
A = {l/k 2 F ){p x ,{ Py ,A}} - {i/4k 2 F )A(d x d y <p), where 
A(r, t) — | A| exp(iip(r, r)) is the center-of-mass gap 
function. J D(p(r, r) denotes an integral over smooth 
("spin wave") and singular (vortex) phase fluctuations. 
Amplitude fluctuations are suppressed below T*. 

It is convenient to eliminate the phase tp(r, r) from the 
pairing term (j^) in favor of d^tp terms [fi — (x,y, r)] in 
the fermionic action. In order to avoid dealing with non- 
single valued wavefunctions we employ the singular gauge 
transformation devised in Ref. |12|: 



V> T -> cxp (itp A )ip h ipi -> exp (1^5)^4, (3) 

where + = Here <Pa(b) is the singular part of 
the phase due to A{B) vortex defects: V x ^ifA(B) = 

2irzJ2i c lid( r ~ r f^)i w it n Qi — ±1 denoting the topo- 
logical charge of i-th vortex and its position. The 
labels A and B represent some convenient but other- 
wise arbitrary division of vortex defects (loops or lines 
in ip(r, t)) into two sets. As discussed in Q this trans- 
formation guarantees that the fermionic wavefunctions 
remain single- valued and the effect of branch cuts is in- 
corporated directly into the fermionic part of the action: 

£ = [dr + i{d r <pA)}^ + i>i\d r + i{d T (p B )]i>i + &H'% 
where the transformed Hamiltonian H' is: 
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with D = (Ao/2fc|,)('?r x ^u + ftyftx) and 7r = p + a. 

The transformation (|3j) generates a "Berry gauge po- 
tential" a M = h(dfj,ipA ~ d^ips) which describes half- flux 
Aharonov-Bohm scattering of quasiparticlcs on vortices 
and mimics the effect of branch cuts in quasiparticle- 
vortex dynamics [l2| This is in addition to the 

"Doppler" gauge field = ^(d^ipA + d^ips) which de- 
notes the classical part of the quasiparticle-vortex in- 
teraction. All choices of the sets A and B are equiva- 
lent - different choices represent different singular gauges 
and Vu is invariant under such transformations. To 
symmetrize the partition function with respect to this 
singular gauge we define a generalized transformation 
(|J) as the sum over all possible choices of A and B, 
i.e., over the entire family of singular gauge transforma- 
tions. This is an Ising sum with 2^' members, where 
Ni is the total number of vortex defects in ip(r,r). 
This symmetrization leads to the new partition function 



Z -> Z = /Z>*t fx>^ JVv^fVa^exp[- f dr f d 2 rC] 
in which the half-flux-to-minus- half- flux (Z2) symmetry 
of the singular gauge transformation (^) is manifest: 

C = &[(d T + ia T )a + iv T a 3 }^ + ¥H^ + C [v^ aj, 
where Co is the 'Jacobian' of the transformation given by 
e -fdrJ d\c I Vip(r, r) (4) 

A.B 

X(5[w M - ^(d^ipA + d^ B )]5[a^ - ^(d^ipA - d^ips)}- 

Here are the Pauli matrices and H = H' . We call 

the quasiparticles & = (V>t, t/'j.) appearing in (^) "topo- 
logical fermions" (TF's). TF's are the natural fermionic 
excitations of the pscudogapped normal state. They are 
electrically neutral and are related to the original quasi- 
particles by the inversion of transformation (pf) . 

To proceed we must extract the low energy, long dis- 
tance properties of the Jacobian (Q). This is done by 
focusing on the fluctuations of two gauge fields and 
a M in the fluid of vortex excitations. We use the saddle- 
point approximation to compute the leading (quadratic) 
terms in £0 for two cases of interest: i) the thermal 
vortex-antivortex fluctuations in 2D layers and ii) the 
space-time vortex loop excitations relevant for low tem- 
peratures (T <C T*) in the underdoped regime (but still 
above the shaded region in Fig. 1). The computation is 
straightforward but the algebra is laborious and will be 
presented elsewhere [Q. Here we quote only the final 
results whose form is ultimately dictated by the symme- 
tries of the problem. For the case i) 
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[(V x v) 2 + (V x a) 2 ] , 



(5) 



where ni is the average density of free vortex defects. 
Both v and a have a Maxwellian bare stiffness and are 
massless in the normal state. As one approaches T sc , 
n i ~ £sc 2 ~~* 0) where t; sc (x,T) is the superconducting 
correlation length, and v and a become massive. Sim- 
ilarly, for the case ii), the quantum fluctuations of un- 
bound vortex loops result in Jl4| : 



Co 



1 

2^ 



K T {d x a) 2 T + ^K t {d * a)\ 



(6) 



where K T ,Ki(i = x,y) are functions of x and T: Ki ~ 
£ sc , K T ~ £J C , z being the dynamical exponent. The 
Maxwell form of Co is dictated by symmetry: the bare 
propagators for t> M and a^, f°(q, ico) and 2?„(q, iu), are 
massless in the normal state and massive within a su- 
perconductor. Note that we dropped from (||) - the 
reason for this is made apparent below. 

The physical picture advanced in this Letter rests on 
the following observations: couples to the TF "charge" 
in the same way as the real electromagnetic gauge field. 
Consequently, if we integrate out TF's in ([!]) to obtain 
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the renormalized (or dressed) gauge field propagators 
X>„(q, iu)) and 2? a (q, icj), we find that f^" 1 (q — > 0, = 
0) — ► const., i.e. the Doppler gauge field is massive. 
This is a consequence of the Meissner response of TF's. 
Physically, this means that the integration over the quasi- 
particles leads to the familiar long range interactions be- 
tween vortices. In contrast, the "Berry" gauge field a M 
couples to the TF spin. This implies that any contri- 
bution of TF's to the stiffness of a M must be massless: 
a singlet superconductor retains the global SU(2) spin 
symmetry ensuring that I?~ 1 (q = 0, iu) = 0) = 0. 

When we combine this with the bare propagators im- 
plied by Eqs. (||||) the following physical picture emerges. 
In the superconducting state, both v M and a M are mas- 
sive by virtue of vortex excitations being bound in finite 
loops. The massive character of i> M and protects the 
coherent TF excitations from being smeared by vortex 
fluctuations. The coupling of TF's to the gauge fields w M 
and a M is irrelevant. This is our inverted Fermi liquid 
phase. 

In the normal (pseudogap) state, the situation changes 
dramatically. The bare propagators for and are 
now massless but the renormalization by the medium of 
TF's screens these bare propagators and still keeps tj m 
massive. Thus, TF coupling to the Doppler shift and 
"spin-waves" remains irrelevant even in the normal state. 
The Berry gauge field a M , however, is now truly mass- 
less since the spin polarization in the medium of TF's 
cannot fully screen the massless bare propagator. In- 
stead, by computing the TF polarization bubble, we find 
V- 1 oc \\/uj 2 +q 2 for (q, iu)) -► |TJ; stiffer than the 
Maxwellian form (||||), but still massless. The massless 
gauge field produces strong scattering at low energies 
and affects qualitatively the spectral properties of TF's. 

The low energy quasiparticles are located at the four 
nodal points of the d xy gap function: (±kp,0) and 
(0, ±kp), hereafter denoted as (1,1) and (2,2) respec- 
tively. Linearizing the fermionic spectrum in the prox- 
imity of these nodes leads to the effective Lagrangian 

C D = Yl *a[ D r ~ ™ f D x (t 3 - iv A D v ax]9 a (7) 

a=l,l 

+ ^ ^aiDr - iV P Dy(J 3 - iv^D x ux]^ a + £o[<v], 

a=2,2 

where ^ is a two-component nodal spinor, a is a node 
index, = + ia^, and Cq is given by (§). We have 
dropped the Doppler gauge field since it is massive 
both below and above T sc and irrelevant for our purposes. 

The Lagrangian ^ and the physics it embodies are our 
main results. In the normal state a M becomes massless 
and the problem of quasiparticle interactions with vor- 
tex fluctuations takes the form of topological fermions 
interacting with massless "berryons", i.e. quanta of the 
Berry gauge field a M . We recognize the above theory as 
equivalent (apart from the intrinsic anisotropy) to the 
Euclidean quantum electrodynamics of massless Dirac 
fermions in (2+l)-dimensions (QED3). The remarkable 



feature of QED3 is that it naturally generates a non- 
Fermi liquid phenomenology for its fermionic excitations. 
This property of QED3 has led to previous suggestions 
that it should be in some form relevant to cuprate super- 
conductivity JTg| , |l7j] . However, the physical content of 
QED3 as an effective low energy theory in this Letter is 
entirely different from those earlier works. 

We now discuss the low energy phenomenology gov- 
erned by the TF propagator: G~ l (k,u>) = G^o(k,a;) — 
S a (k, uj), where G~q is a free Dirac propagator at node 
a. We first consider the T = case in the isotropic limit 
(vp = va) where explicit results are readily obtained. In 
vpk x a 3 — v&kyO~x and we find 



this case G q q = uj 



3tt 2 N 



{-ui + Vpk x a 3 + v^kyai) ln(A/p), (8) 



with p = (- 



'Jpk x 



+ V 2 A 



h 2 ,) 1 / 2 , A a high energy cutoff, 



and N = 2 is the number of pairs of nodes. 

The essential feature of the TF propagator is the sin- 
gular behavior of the self-energy £ Q (k, u>) which arises 
from the massless nature of the dressed berryon propaga- 
tor 2? a (q, w) and is logarithmic in the leading order. This 
result can be formalized as the leading term in a large N 
expansion. Ultimately, the resummation of such an ex- 
pansion ]l6[ | yields a power law singularity G a cx j/' -1 , 
with a small exponent 77 = — 8/3ir 2 N [[TJ|. For our pur- 
poses, having to deal with both the anisotropy and finite 
T, the leading order form (0) is more convenient since 
it allows for explicit computation of various quantities. 
Once we move beyond the leading order, the vertex cor- 
rections to E Q are necessary and the algebra becomes im- 
penetrable. Furthermore, the available experiments are 
unlikely to distinguish between 77 = + and a small finite 
exponent. 

The singularity in E Q heralds the breakdown of the 
Fermi liquid behavior in the normal state. To see this 



consider Eq. (Bh for E k = 



find E Q oc -(8/37r 2 A^)wln(A/V^ 
fermion pole vanishes as ui 
its width goes as SJJ 



lk 2 v fl 2 « \u\. We 

. The residue of the 
^ 1/ In while 
This behavior is 



-uj 
0, Z(w) 
-(4/3ttA0M. 

reminiscent of the MFL expression for the self-energy, 
assumed on phenomenological grounds by Varma, Abra- 
hams and collaborators Pq| . Note, however, that our 
QED3 TF propagator implied by Eq. (^) remains qual- 
itatively different from the MFL Ansatz 15 , both by 



the fact that ln(A/p) is replaced by a weak power law 
(thus algebraic Fermi liquid) and by the momentum de- 
pendence of S Q (k, u>). As shown below it is this com- 
bined momentum-frequency dependence that provides a 
natural explanation for some of the remarkable features 
of the fermionic spectral function in cuprates observed in 
the ARPES experiments JT9|. Also, we emphasize that 
our results apply to the pseudogap phase below T*. The 
physics of the normal state at higher temperatures is be- 
yond the scope of our present theory. 

Inspection of Eq. (Q) reveals that E Q has imaginary 



part only inside the cone defined by uj > v F k 
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FIG. 2: Energy versus momentum distribution curves of 
A(k, oj). Left: EDC cut taken for k = (coincident with 
a nodal point), and MDC cut taken for u> — 0~ and k y — 
0. Both curves have been broadened (by the same amount) 
to simulate the finite resolution of an ARPES experiment. 
Right: The corresponding spectral function density plots for 
isotropic (top) and anisotropic vf/va = 17 (bottom) cases; 
to be compared with Figs. 1 and 2 of Ref. p9(] . 



outside this cone vanishes. This implies that TF spec- 
tral function plotted as a function of momentum at fixed 
u> (MDC) will be very sharp close to the Fermi surface, 
while the corresponding energy distribution curve (EDC) 
will be broad. This is illustrated in Fig. || where we plot 



the spectral function A(k,u) 



L Im[G Q (k, W )]n de- 



duced from Eq. (|S|) . We note that precisely such striking 
asymmetry between the EDC and MDC cuts is observed 
in the ARPES data |l9). 

These qualitative features of the spectral function sur- 
vive at finite temperature and away from the isotropic 
limit. Unfortunately, away from this simple limit the 
precise form of the TF propagator is not known: as 
soon as the "relativistic" invariance of the T = prob- 
lem (^) is lost, analytic calculations become intractable. 
We find that for T <C oj, the self energy retains its 
T = form Eq. (^J) with a small temperature correction. 
On the other hand when T > to,E k we find — T, 
qualitatively consistent with the original MFL conjecture 
£^ ~ max(cL>,T). We note that such T-linear scattering 
rate has been deduced from ARPES experiments [|l9| . 

In ARPES one measures the spectral function of 
real electrons not of TF's. While the inversion of the 



transformation (JsJ) after the phases have been coarse- 
grained and replaced by the gauge fields is a daunting 
task, our theory insures the gauge invariance with re- 
spect to <x M of the true electron propagator. The sim- 
plest such gauge invariant propagator is Gf 1 1 oc (x, x') fts 
(exp(i/* ds M o At )[#(a;)#t( a ;')] u ) ) where x = (r,r). By 
employing a gauge in which the line integral of van- 
ishes [^0j , we have computed the asymptotic behavior of 
G clcc (x,x'). We find Q that it exhibits a power law 
singularity with the exponent r( = 2n — — l6/3ir 2 N. 
This strongly suggests that the true electron propagator, 
whose precise form within QED3 is unknown at present, 
will exhibit a power law with small positive exponent. 

In conclusion, we argued that the pseudogap regime 
in cuprates can be modeled as a phase disordered d- 
wave superconductor. Such assumption naturally leads 
to a QED3 theory for the massless Dirac "topological" 
fermions interacting with a massless gauge field of vor- 
tex "berryons" . Coupling to the massles gauge field de- 
stroys the Fermi liquid pole in the fermion propagator 
and generates algebraic Fermi liquid. Lacking any en- 
ergy or length scale this theory can be thought of as being 
critical independently of the actual doping level x. Below 
T* the low energy spectral properties of the fermions are 
therefore regulated by a quantum critical line. In this 
regime the low energy fermiology, including thermody- 
namics, transport, and density and current responses are 
all controlled by the universal properties of topological 
fermions and vortex berryons encoded in the anisotropic 
QED3 Lagrangian (Q) . Eventually, this peculiar quantum 
critical behavior gives way to the actual superconduct- 
ing phase at T sc (x) and the Fermi liquid character of the 
nodal quasiparticles is restored as vortices bind into finite 
loops. At very low doping, hole Wigner crystal, SDW, 
and other low T phases become possible, reflecting the 
strong Mott-Hubbard correlations. 

The authors are indebted to G. Baskaran, M. P. A. 
Fisher, S. Girvin, I. Hcrbut, C. Kallin, S. A. Kivelson, A. 
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blay, O. Vafek, P. B. Wiegmann and A. Zee for helpful 
discussions. This research was supported in part by NSF 
Grants PHY99-07949 (ITP) and DMR00-94981 (ZT). 
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